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ABSTRACT 

HQET currents with the quantum numbers of the ground-state 
baryons are discussed. One-loop anomalous dimensions are cal- 
culated, and exact one-loop matching to QCD currents is found. 
Two-point correlators of these currents are calculated taking into 
account d < 9 terms of the OPE. Sum rules for heavy baryons 
Aq and Eg are analyzed. Three-point correlators of two bary- 
onic currents and a heavy-heavy velocity changing current are cal- 
culated with the same accuracy. The baryonic Isgur-Wise form 
factors are estimated from the corresponding sum rules. 
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1. Introduction 



It was noticed long ago Q that properties of hadrons with a single heavy quark 
are simple: it's mass (and hence flavor) and spin orientation are irrelevant 
to the leading order in 1/to. Recently this qualitative physical picture was 
incorporated into the formal framework of the Heavy Quark effective Theory 
(HQET) |2,p. The heavy quark spin symmetry Q leads to relations among 
meson [Q, 5 and baryon [6-8] form factors. An elegant general discussion 
can be found in . Not only the orientation but also the magnitude of the 
heavy quark spin is irrelevant, so this symmetry is extended to the superflavor 
symmetry 10 . 

Nonperturbative methods (such as sum rules [pd]| ) are needed to obtain 
hadron properties in HQET. HQET sum rules for mesons and baryons were 
first considered in jl] . They should be improved in order to take into account 
the currents' renormalization in HQET. One-loop anomalous dimensions of 
HQET meson currents were found in , and two-loop ones — in . HQET 
meson sum rules with the proper account of these effects were considered 
in . the meson Isgur-Wise form factors were considered in the frame- 
work of three-point sum rules . 

Baryonic currents and sum rules in QCD were considered in a number of 
papers [17-20]. One-loop anomalous dimensions of QCD baryonic currents 
were found in [ pT[ , and two- loop ones — in [ |2^ (where the perturbative cor- 
rection to the sum rules was also obtained) . Three-point sum rules for baryon 
form factors were discussed in [^3| . 

In the present work we discuss the HQET currents with the quantum num- 
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bers of ground-state baryons. In Sec. ^we calculate their one-loop anomalous 
dimensions (as was done for mesons in jl^), and obtain the exact one- loop 
matching of these HQET currents to QCD ones (as was done for mesons 
in 1^). Note that non- logarithmic terms in the one-loop matching are useless 
unless the two-loop anomalous dimensions are known. They can be found 
using the methods of [ p^ ; we hope to return to this question later. 

Our further analysis of the heavy baryons is similar to the detailed light 
baryon analysis Q. In Sec. || we consider two-point correlators of the HQET 
baryonic currents. OPE for diagonal correlators contains even-dimensional 
terms; perturbative terms (d — 0), (G^) corrections {d = 4), (qq)'^ terms 
{d — 6), and corrections to it (d = 8 . . .) are taken into account. OPE for 
nondiagonal correlators contains odd-dimensional terms starting from (qq) 
{d — 3); ttiq and (G^) corrections {d — 5 and 7), and {qq) terms {d — 
9) are also taken into account. After that, the sum rules for Aq and Eg'' 
following from these correlators are analyzed. The diagonal sum rules have 
large continuum contributions because the spectral density grows like u^. In 
the nondiagonal sum rules it grows like , and continuum contributions are 
not so important. But higher power corrections d > 7) are poorly known. 
Nevertheless, both types of sum rules give good results consistent with each 
other. We disagree with some results of 

In Sec. ^ we consider three-point correlators of two HQET baryonic cur- 
rents and a heavy-heavy velocity changing current. We calculate diagonal 
and nondiagonal correlators with the same accuracy as in the two-point case, 
and obtain the baryon Isgur-Wise form factors from the sum rules. 



2. Baryonic currents in HQET 

Hadrons in HQET are classified according to the light fields' angular momen- 
tum and parity j'^ . For the ground-state baryons, light quark spins can add 
giving = 0"*^ or 1'*'. In the first case their spin wave function is antisymmet- 
ric; Fermi statistics and antisymmetry in color require antisymmetric flavor 
wave functions. Hence light quarks must be different; if they are u, d then 
their isospin / = 0. In the second case the flavor wave function is symmetric; 
if light quarks are m, d then their isospin 1=1. This gives us the baryon 

with 7 = called Ag, and the degenerate and |^ baryons with 1=1 
called Eq and Eg. 

Baryonic currents have the form J = e"'^'^ {qf'^ CTTq2)T' Q'^ , where C is 
the charge conjugation matrix, means q transposed, r is a flavor matrix 
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(symmetric for Sg and antisymmetric for Ag), and Q is the effective field 
satisfying 70Q — Q. We sliall abbreviate it to J = {q'^Crq)r'Q. A light 
quark pair with j'^ = 0+ corresponds to the current a = q^C^5q, and with 
j"^ = 1+ — to a = q^C^q. One can easily check it using the P-conjugation 
9 ^ 709- It is also possible to insert 70 into these currents without changing 
their quantum numbers. The current J= aQ has spin 1/2; the current J= aQ 
contains spin 1/2 and spin 3/2 components. The part J3/2 = .7+ ^77 • J 
satisfies the condition 7 • J = and hence has spin 3/2. The other part 

Ji/2 -577 • J== I775J1/2, Ji/2 = a • 775Q has spin 1/2. Finally we arrive 
at the currents 

Jai = {q^Cj5q)Q, JA2 = (g^C7o75(?)Q, 

Jsi = (g^Cfq) • 775Q, Js2 = (g^C7o7g) • 775(9, (1) 

Je-i = ('?'^C7g)Q + 17(9^^79) • 7Q, 

Js'2 = {q^Cjojq)Q + \l{q^C-io^q) ■ 'jQ. 

This classification follows ideas of the currents for Ag and Eg first 

appeared in 0|. 




a b c 



Figure 1: One-loop proper vertex 

Now we consider the one-loop renormalization of these HQET currents. In 
the MS scheme (the space dimension D ~ 4 — 2e) the bare fields are related to 
the renormalized ones by qo — fi^^ Zqq^ Qq — Zq Q, where pr = ^^jf-, 
/Lt is the normalization point. In the Feynman gauge the renormalization con- 
stants are Zq = Zq = l-h2CF^, where Cp = TV = 3 is the 
number of colors. The bare current Jo = qoqoQo = ZjJ, or qqQ = ZrJ, 
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where Zj = ZqZg^ Zr- Then the matrix element F ~ (o\qqQ\qqQj = 

Zr(o\j\qqQ^^ where the matrix element of J is finite. Calculating the ul- 
traviolet divergence of the vertex T Fig. ^ in the Feynman gauge, we obtain 
Zr = 1 + (4^ + 2) , where j^^Tj^ = HT, Cb = Using and 

Zq we find Zj = Zp; hence the anomalous dimensions = '^^"^^^l are 

1.--Cb^(^ + 2\, (2) 



' 27r V 4 

„ as ~ ~ i as ^ ^ "s 

7A1 = — otyB , 7A2 = 7S1 = -TT^B , 7S2 = — C-B . 

TT Z TT TT 

They depend only on the light quark part of the current. 

Now we are going to establish relation between these HQET currents and 
QCD ones. A QCD current {q^CTq)T'Q matches the corresponding effective 
current Aj if they give the same physical matrix elements. In order to cal- 
culate on-shell matrix elements, we have to use the on-shell renormalization 
scheme in which propagators in the on-shell limit are free. For the "massless" 
fields g, Q, the bare on-shcU propagators get no corrections because loop in- 
tegrals are no-scale (ultraviolet and infrared divergences cancel). Therefore 

— 1/2 

the on-shell renormalized fields coincide with the bare ones: q — Zq qos, 
— — — 1/2 — 

Q — Zq Qos- Note however that although the expressions for renormaliza- 
tion constants Zq, Zq are the same, all divergences in them are infrared ones 
because these Z factors relate renormalized (ultraviolet-finite) fields. For the 

massive quark field Q we have Q = Zq^^^Qos, Zq = 1 + Cpj^ (| — 3L + 4), 
2 

where L = log Note that the infrared divergence of the on-shell massive 
quark propagator Zq is the same as that of the effective quark propagator 

Zq- 

We have the currents j = Z^^ Z;^^^ Z^^^^ qosqosQos, J = Z^^Z^'^Z^^^^ 
9os9osQos- Hence the on-shell matrix elements are {0\j\qqQ) — 
Z^^Zq-^Zg^^^T, (0\j\qqQ) = Zf^Z-iZ^^^^F, where the proper vertices F, F 
are depicted on Fig. 0. Hence we obtain the matching constant 

1/2 

r/^r . Zq\ 1 , ^ / 3 




A^Aq..^, Aq^\-^\ =i + Cf—[-L 
^F/Zr' [Zq "^4^12 



■^of course, it is a pure formality to write N in formulae based on the fact that a baryon 
contains three quarks. 
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Here ultraviolet divergences cancel in T/Zy, T/Z-p by definition; infrared di- 
vergences cancel between these two expressions because the infrared behavior 
of QCD and HQET is identical; Aq is finite. We choose all quark momenta 
in HQET to be zero; this corresponds to the heavy quark momentum mv 
{v = (1,0)) in QCD. Then the HQET loops (Fig. |) vanish. The diagram 
Fig. ^ vanishes in F (but not in Zr). Logarithmic terms in the matching con- 
ditions are determined by the difference between QCD and HQET anomalous 
dimensions. 

The choice of QCD currents corresponding to the HQET currents (|l|) 
is not unique. We restrict ourselves here to the simplest variant with F' — > 
If we denote {{qTc[a^,,T]q)^^T'i^a^,Q) = f{D) (j), then 

[ + Csfj (|(L - 3) + i^)] . Indices v 

are purely space-like, hence the matching conditions are the same for pairs of 
currents with F — > 70F. Finally we obtain (in the scheme where 75 anticom- 
mutes with all 7^) 



1+70 p/ 1+70 
2 2 



the matching constant A — A, 



1 



-70 



-('Z^C7^9) 
1 



1 



- 70 1 
5 — 7^75 — 



-70 



Q^Aq 



Cb^{2L~4) %u (3) 



D - 1 



-70 



Q 



These equations should be used to obtain the QCD matrix elements from the 
HQET ones calculated, for example, from the HQET sum rules (Sect. 0). 



3. Two-point correlators 

Two-point correlators have the general structure (J= {qTr^ C^^lf' Qt' ) 

^ (Tji{x)%iO)) = (^F'ii±^r;) S{x)2TrTT+Tnixo). (4) 

The expression (||) without the first factor is the correlator in which the 
heavy quark spin is switched off (its propagator is S{xq) — —i'd{xo))- For 
both — 0+ and 1+ there are two diagonal correlators and one nondiagonal 
one (see (|l|)). Correlators for the physical spin symmetry multiplets (Ag, En, 
Eg) are expressed via these ones (||) due to the superflavor symmetry [lOfl . 
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The normalizing factor is T — ■iTrrir2 for diagonal correlators and T = 
|;Tr7orir2 for nondiagonal ones. 

OPE for diagonal correlators contains even-dimensional terms. We take 
into account the perturbative {d = 0) term (Fig. ||a), the gluon condensate 
(d = 4) term (Fig. |b), the {qqf (d = 6, 8. . . ) terms (Fig. §;). The diagram 
Fig. ||d gives a, d = 6 contribution; it is of the order of unknown perturbative 
corrections to the diagram Fig. ||c, and is not taken into account. We use the 
fixed-point gauge Xp,Afj_{x) = in which the heavy quark does not interact 
with gluons. The methods of calculation of correlators in this gauge are 
reviewed in l24|. We obtain 



Nl S{t) 



1 + c 



32iV(A^- 1) 



(5) 



The terms Ila+b appear to be the same for jiji and J2J2 correlators; c = 1 for 



An and 



We have factorized a four-quark condensate in He into 



i-Q ana — i for Sg''. 
two two-quark ones. In this approximation He is also the same in two diagonal 
correlators (though the factorization approximation is thoroughly checked 
only for products of vector and axial currents). These correlators can't strictly 
coincide at least because Ji and J2 have different anomalous dimensions (||). 

The nonlocal quark condensate {q{0)q{x)) = (qq) ^1 + ^-^^ — \- ^"9^^ ^ + 
where in the last term the factorization is used for d = 7 quark-gluon 
condensates. The Gaussian anzatz {q{0)q{x)) = (qq) exp ^ ™°g ^ 



was pro- 

^ term in it is about 3 times larger than in the 



posed in |25[] instead. The x 
factorization estimate. 

The perturbative {d = 0) and (qq)'^ {d — 6) terms were first considered 
in [0. We agree with the perturbative result. The {qq)'^ term in jl] had 
different signs for An and Eg while these signs are the same in (^). Because 
of this, the results ||] on Eg and its difference from Ag are incorrect. 

OPE for nondiagonal correlators contains odd-dimensional terms. The 
diagram Fig. Bd gives d = 3, 5, 7. . . terms; Fig. ||f — d = 7. . . ones. The 
diagram Fig. He contributes a d = 9 (qq)'^ term. We obtain 



Nl {qq)Sit) 



16 



96iV 
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^! (qq) Sjt) 
16n'^N{N- l)t 



6N 



(6) 



.CpNlnas {qq)S{t)t^ 



144iV3 



The factorization approximation is used for d = 7 quark-gluon condensates, 
so these terms are order-of-magnitude estimate only. In this approximation 
the diagram Fig. ||f does not contribute. 

Correlators obey the dispersion representation 



n(c.) 



p{e)de 



e — Lo — 



- + •••, n(t) = -S{t) / p{uj)e-'^'dLo + • • • (7) 



A subtraction polynomial in Ii-{uj) (denoted by dots) leads to 5{t) and its 
derivatives in n(t). We analytically continue correlators from t > to imag- 
inary t = —IT. Then H{t) and p{uj) are related by the Laplace transform 



a+ioo 



n(r) = I J p{uj)e-'^^dLO, p{u) = ^ J n(r)e"-rfr, 



(8) 



where a is to the right from all singularities of n(r). A term S{t)/t'^ [n > 1) 
in U{t) gives i"d{uj)ui"~^ /{n— 1)! in p{u;); a term S{t)t" {n > 0) corresponds 
to (— (w). We obtain for the diagonal and nondiagonal correlators 



5! 



'■i2N{N ~ 1) 



Pn{uj) = 



N\ {qqf 
47V2 

Nl {qq 



Cp-n-^cts {qqf 

144iV2 



5"'{u) 



(9) 



6iV 



We have also verified the leading terms of these formulae by a direct appli- 
cation of the Cutkosky rules in the momentum space. 
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We equate the OPE-based expressions for II{t) to the dispersion 

representation (|^) with a model spectral density containing the lowest baryon 
contribution and the continuum contribution. This procedure is equivalent 
to the nonrelativistic Borel sum rules 

Baryon residues in QCD are usually defined as {0\j\B) = fus where the 
relativistic normalization of the baryon state \B) and its wave function ub 
is implied (and Tr r^r = 1 is assumed). This normalization is senseless in 
HQET; one should use the nonrelativistic normalization instead. Both sides 
of this equation are divided by ^/2m, and its form does not change. So 
baryonic residues / are constant in the heavy quark limit up to logarithmic 
renormalization effects. We have the lowest baryon's contribution to the 
spectral density p{w) = ^S{uj — e) where e is the HQET baryon's energy 
(i. e. mass minus the heavy quark mass). 

We adopt the standard continuum model: its spectral density is equal to 
the theoretical one starting at a continuum threshold ujc (of course, S^"''{uj) 
terms don't contribute). The continuum contribution is thus subtracted from 
the theoretical part of the sum rule leaving an integral over the resonance's 
duality interval (up to uJc)- A term S{t)/t^ after subtracting the continuum 
contribution gets the factor fn-i{uicT) where /n(x) = 1 — e-~^Y^^=oln\- 
Therefore the formulae (^ are not necessary to construct the sum rules. We 
include all S{t)/t^ terms with n > 1 to the spectral density, as was done 
in ||l^. If power corrections are not included, it means the absence of such 
subtraction for them what leads to their overestimation. 

The two diagonal correlators give the identical sum rules for (0|Ji|S) and 
(0|j2|i?). Therefore these matrix elements are equal with the accepted ac- 
curacy. Following ||l[, we introduce the dimensionless variables ''' — 

P = ^^^n, mo = 4kEo, ,Z\^f_\^ = {kEG)\ = fci?,, e = fcS,, 
— —jf^ (qq)- Then we arrive at the two sum rules 



exp ( -2^ ) (10) 



= 2E^ 



M ) - (1 - 2) ;^/o ( ) + 3 (1 - 2) ^ + 



We have used the Gaussian anzatz for the nonlocal quark condensate in the 
diagonal sum rule. The nondiagonal one is used at larger E where the nonlo- 
cality is not so important. For consistency, we use the factorization estimate 
for all of the diagrams Fig. ||d-f (including the nonlocality) . 

Note that (0|Ji|i3) and {0\j2\B) can't be always equal because they have 
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different anomalous dimensions; hence the anomalous dimensions effects are 
out of our control here. We obtain these matrix elements normalized at 
a typical Borel parameter scale; if one wants to have them at a different 
normalization point (e. g. at mg for matching with QCD), one should use 
the anomalous dimensions (|2|). 

At the standard value of the quark condensate, the energy unit k = 
280MeV. It is stated in ||l8| that light baryons are described better if the 
quark condensate is reduced by 20%. This leads to fc = 260MeV; this value 
was also used in Q. We also accept it; the standard values of ttiq and (G^) 
then give Eq = 0.85, Eq = 0.60. 

The results of the diagonal sum rules analyses are shown in Fig. ^. In the 
selected range of the Borel parameter E, the perturbative contribution, the 
gluon condensate one, and the quark condensate one constitute for Ag 60- 

40%, 25-15%, and 30-35%, correspondingly. For Sg\ the gluon condensate 
contribution is 3 times smaller and has the opposite sign. The continuum 
contribution (which is subtracted from the theoretical side of the sum rule) 
rapidly grows, and is several times larger than the total result at the right 
end of the interval. This means that the relative error due to the rough con- 
tinuum model is multiplied by 1-5 in the total result. On the other hand, the 
contribution of poorly known higher power corrections (with d > 8) becomes 
large at the left end of the interval. Namely, the nonlocal quark condensate 
is several times smaller than the local one at this point, and the d — 6 con- 
tribution is almost completely compensated by d > 8 ones. The remarkable 
stability of the sum rule in the lower part of this interval is a plausible argu- 
ment in favour of the Gaussian anzatz for the nonlocal quark condensate (the 
results with, i. g., the factorization anzatz are not so stable in this region). 
Thus we obtain 

Aq : e = 780MeV, lu^ = 1200MeV, / = (1.8 — 2.7) • IQ-^GeV^, 
S^*' : e = 990MeV, lu^ = 1460MeV, / = (2.9 4.1) • lO^^GeV^ 

Taking into account the recent experimental result |2^ ttia^ = 5640 ± 50 ± 
30MeV, we obtain a reasonable value mi, = 4860MeV. The Af,-Eb splitting 
of 210MeV is in agreement with the potential model expectations of 190MeV 
(see e. g. [^^). The distances from the resonance energies to the continuum 
thresholds in Af, and Sft channels are 420MeV and 470MeV; they approxi- 
mately correspond to the distances to the first excited states in these channels 
expected in the potential model ^ 460MeV and 405MeV. 

The results of §] for Aq are s = 700MeV, / = 2 • lO^^GeV^; it is an 
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2 



1 




Figure 3: The diagonal sum rules, a) Er as a function of E at various Ec from 
the logarithmic derivative of the sum rule: three lower curves — results for Aq 
at Ec = 4.1, 4.6, 5.1; three upper curves — results for Eg at Ec = 5.1, 5.6, 
6.1. b) n as a function of E at various Ec, Er: three lower curves — results for 
Aq at Ec = 4.1, Er = 2.75: Ec = 4.6, Er = 2.95; Ec = 5.1, Er = 3.2; three 
upper curves — results for Sg^ at Ec — 5.1, Er — 3.6; Ec = 5.6, Er = 3.8; 
Ec = 6.1, Er = 4.05. 
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order of magnitude too lowQ. The sum rule for Sq in [Q is incorrect. Ac and 
Sc were considered in the second paper of [ pO[ using relativistic sum rules. 
Two matrix elements obtained there reduce to (0|Jai|Ac) in the heavy quark 
limit, and give /a^i = 3.2 • lO^^GeV'^ and 4 • lO^^GeV"^; one matrix element 
reduces to (0|Js2|Sc), and give /s,2 = 7.2 • IQ-^GeV^. 

The results of the nondiagonal sum rules analyses are shown in Fig. ^. 
In the selected range of the Borel parameter E, the mg correction in the 
Aq sum rule reduces starting from 40%; in the case of Sg-* it is 7/3 times 

larger. Therefore the left end of the interval is not very reliable for Sg . The 
estimated contribution of poorly known d > 7 corrections docs not exceed 
10%. The continuum contribution grows not so strongly as in the diagonal 
case, allowing us to use larger values of E. 

The sum rule for Ag prefers somewhat lower values for the effective con- 
tinuum threshold than in the diagonal case, and gives the same results for 
the resonance energy and the residue. The sum rule for Eg gives somewhat 
too high values for the mass and the residue. In general, the agreement of 
these two completely independent sets of sum rules gives us more confidence 
in the results. 



4. Three-point correlators 

Now we shall consider correlators of two HQET baryonic currents Ji,2 — 
{q'^CT I ^2C[)^'i 2Q1.2 (where Q1.2 = i^i^2Qi,2 are the effective heavy quark fields 

with the velocities ui,2, vi-V2 = ch (p), and a heavy-heavy current J ~ QiTQ2. 
They have the structure 

{TMx,).m%ix2)) = (r;^^ri±^r;) 2Tr rr+ 

00 00 

J dti5[xi ~ Viti) j dt25{x2+V2t2)K{ti,t2). (11) 


The expression (^l|) without the first factor is the correlator in which the 
heavy quark spin is switched off. For j"^ = 0+ K{ti,t2) is a scalar function; 
for = 1+ it is a tensor Kf^^, = K\\ei\\fj_e2\\„ + K^5±^i,, where ei|| = (w2 - 
chifvi)/ ship, 6211 = —{vi — ch (y5W2 ) / sh iy9 are the light fields' polarization 

*At the quoted value of /aq , the right-hand side of the equation (45) in jlj is about 
two orders of magnitude less than the right-hand side. 
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10: 

0.9 1 1.1 1.2 1.3 1.4 1.5 

70r 
60 : 

50 ; 

40 ; 

30 ; 
20 ; 

10; 

1.75 2 2.25 2.5 2.75 3 3.25 3.5 

Figure 4: The nondiagonal sum rules, a) Er as a function of E at various Ec 
from the logarithmic derivative of the sum rule: two lower curves — results for 
Aq at Ec = 3.4, 3.0; three upper curves — results for Eg^ at Ec = 5.1, 5.6, 
6.1. b) n as a function of E at various Ec, Er: two lower curves results for 
Aq at Ec = 3.4, Er = 2.65; Ec = 3.9, Er = 2.95; three upper curves — results 
for S[!^ at Ec = 5.1, Er = 4.2; Ec = 5.6, Er = 4.4; Ec = 6.1, Er = 4.65. 
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vectors in the scattering plane, and S^^i, = ^ ei j_^e2±,y = [ch ifi{vi^V2v + 
V2piViu) — vi^viu — V2nV2v\/ sh^ </? — g^i,, eii_ — e2± are two orthonormalized 
polarization vectors orthogonal to this plane. In each of the three cases (0+, 
there are two diagonal correlators and one nondiagonal one. 
In the limiting case (p — we have 



K{tut2)^n{ti+t2) 



(12) 



Indeed, let's consider any diagram for the two-point correlator in the coordi- 
nate space (for simplicity, with the scalar heavy quark). Vertices along the 
heavy quark line have times ta < ti < ■ ■ ■ < i„_i < and integration 
in ti,. . . , tn is performed. Consider the diagrams obtained by inserting the 
heavy-heavy vertex (with the time t and (p = 0) to all possible places. These 
diagrams have the integration regions to < ti < ■ ■ ■ < tm-i < t < tm < 
• • • < tn-i < tn- These regions span the whole integration region of the 
original diagram, therefore the sum of all these three-point diagrams is equal 
to the two-point diagram. In the momentum space the equation (12) reads 



K{uJi,UJ2) 



For Ttq correlators at (yS = we have K\\ = K± 



because all directions orthogonal to v are equivalent 
For the diagonal correlators we obtain (Fig. |^) 



Ka+b{ti,t2) ^ -i 



N\S{h)S{t2)a 



7r4(a;2)4 

^^^i^(4(A^-2)^?t^hV + 3c(.Y) 
N\S{h)S{t2) ,_, 



Kc{ti,t2)-i 



(13) 



47V2 



{qi-V2t2)qiviti)y 



where a; = witi -1-^2^2, = tl + t2 + 2chiptit2; a 



for Al, S||2, S^l cases 
^ ~ chip{tl+tl) + 2tit2 for Aa, Eiil, 



(1, 2 refer to ji^2 correlators), and a = y 

S^2 cases; c = 1 for Al, c = i {l + 2ch(/j|^^ for A2, c^\[l - 2ch(^|^^ for 

I]||l, c -i for S||2, E_l1, S±2. These results agree with (|) at ^ = (0). 

The nonlocal quark condensate {q{x)q{y)) in the 0-gauge (i. e. the 
fixed-point gauge (a; — xo)^A^{x) = with xq — 0) can be rewrit- 
ten in the gauge-invariant form {q{x)E{x,0)E{0,y)q{y)) , where E[x,y) — 
P exp{—i) Afi{z)dz^. Using the translational invariance we rewrite 
it as {q{0)E{0,-x)E{-x,y - x)q{y - x)), or {q{0)E{-x,y - x)q{y - x)) 
in the new 0-gauge. In the factorization approximation (giving 
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only a rough estimate) we can replace E(x, y) in this formula by 
{j^TTE{x,y)) = 1 + !^^i^((xy)2 -xV). This gives {q{-V2t2)q{viti)) = 

(g(O)g(x)) 1 + ^-^^ ^ 

For the nondiagonal correlators we obtain (Fig. |^) 
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a 1 



16 



+ b 



tit2 sh ip 
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{{x^f + 2t\tl sh^) 



(14) 



Keitl,t2) 



NlSih)S{t2){qq) 
16tt^N{N - 1) 
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Kc{h,t2) 



CFNlS{ti)S{t2)TTas {qqy {h + chipt2)d 



where a ^ ti + ch(pt2, b — t2, d ~ I for A and a = t2 + chipti, b = ti, 
c? = ch for cr|| ; c = 1 for A. c = — i for S^. These results agree with (^ 
at V3 = (|l|). 

Three-point correlators obey the double dispersion representation 



K{uJi,UJ2) 



p{£i,£2)dside2 



[Si - UJi - iQ){£2 - UJ2 - 

K{hM) = S{h)S{t2) f p{L^i,L02)e-'^'''-"^'''dujidLj2 + 



(15) 



Subtraction terms in K(uji, LU2) (denoted by dots) are polynomials in uji with 
coefBcicnts that are arbitrary functions of UJ2 (given by single dispersion inte- 
grals), or vice versa. These terms give 5^"^(ti) times arbitrary functions of t2 
(or vice versa) in K{ti,t2- We analytically continue K{ti,t2) from ti^2 > to 
ti_2 = — iTi_2- then K{ti,T2) and p{ui,uj2) are related by the double Laplace 
transform 



K{ti,T2) = - J p{LUl,0J2)e-^'' 



^dwidw2, 



(16) 
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a+ioo a+zoo 



a — too a— ioc 

All considered correlators have the form K{ti,T2) = P(ti,T2)/(— x^)", where 
P{ti,T2) is a polynomial, and —x'^ = + r| + 2chipTiT2. It is convenient 
to introduce the new variables zi = e'^/^Ti + e^'^/^r2, Z2 = 6*^/^x2 + e^'^'/^ri, 
r^i = e'^/^wi - e-™'^P''V2^2, ^2 = ev/^LU2 - e-'^'^'P'^'/'^LUi. We have -x^ = 
Z1Z2, UJ1T1+UJ2T2 = {i^iZi + fl2Z2)/2sh(p, and integrals for p{u}i,uj2) factorize. 
They don't vanish only in the region Sli > 0, fl2 > 0. This region has the 
form of a wedge e~'^ < ^ £ e*^ (Fig- §)■ Therefore it is convenient to use 
the parameterization UJ12 — Lu{l±r]t]i ^) (or f^i 2 = 2ix' sh ^(1 ± 77)) in which 
the physical region is —1 < 77 < 1. It is evident from ( [l5| ) that in the limit 
(fi ^ K(ti,t2) depends only on ti + 12 as it should (p^), and is given by the 
single dispersion representation (^ with 

+1 

p{lo)= mii2tjthf j p(u;(l + ??thf),w(l-rythf))d?7. (17) 



Figure 6: Physical region of double spectral densities 
For the diagonal correlators we obtain 



27r4 sh ip 



u)^A{7-i) 7ra,(G2) (N-2 



5!ch*f 32iV(Ar- 1) V6sh> 



-B{^)+C{i^) 



(18) 



where A(r/) = f (1 - 'q'^f, B{t]) = S'{1 + tj) + d'{l - ??) - (3 + Bsh^ f ) 
(<5(l + r/)+5(l-?7)-l) + Bsh-^f for Al, S||2, E^l; A{r]) = |(1 - ??■*), 
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Birj) = i (,5"(1 +Tj) + S"{1 - rj)) - (l + ish^ |) (,5'(1 + 77) + S'{1 - 77)) + 
2 (1 + 2 sh^ I + 4 sh* I) ((5(1 + r;) + 5(1 - 7?) - 1) + 8 sh'' I for A2, S,, 1, S_l2; 
C(?7) = 1 for Al, C{r]) = i ((5(1 + r?) + 5(1 - 77) + 2ch(p) for A2, C{r]) = 
i((5(l + 77) + 5(l-77)-2ch(^) for I]||l, C(77) = -i for I]||2, S^l, C(?7) = 
— i (5(1 + r;) + 5(1 — 77)) for T,±2. For the nondiagonal correlators we obtain 



{qq)uj 



l±77th| (19) 



27r2sh(p 

- (1 - f ) {e^^'Sil T 77) + e--/^5(l ± 77)) 

In this formula we assumed A'^ = 3 because it becomes much simpler in this 
case. Upper signs are for A, E^, lower sign — for E||; c = 1 for A, c— — -i for 
E||_L. Terms with two 5-functions nonvanishing only at the origin uji = 0, 
llJ2 = are omitted in equations (p^[-p^. These formulae agree with (||) at 
if — (p7|). We have also verified the leading terms in these formulae by a 
direct use of the Cutkosky rules in the momentum space. 

In order to obtain the information on the lowest baryons' form factors, we 
equate the OPE-based expressions (p^-p^) for K{ti,T2) to the double disper- 



sion representation (15) with a model double spectral density containing the 
lowest baryons' contribution in both channels, continuum contribution, and 
the mixed contributions with the lowest baryon in the one channel and con- 
tinuum in the other one. The last contribution is exponentially suppressed at 
sufficiently large ri,2, and is neglected as usual Q. The subtraction terms 
don't contribute at n ^ 0, T2 ^ 0. The applicability regions of the sum rules 
are symmetric with respect to the interchange of ri and T2 (or nearly sym- 
metric in the nondiagonal case). Therefore we shall not loose an important 
information if we restrict ourselves to the diagonal ri = T2 = t/2. 

The contribution of the lowest baryons in both channels to K{ti,t2) has 
the form i/i/|e~*^i*i~'^2t2^((,h(^), where ^(ch (ys) is a scalar function for j'^ = 
0+ and — C| 1 6111^162111/ + for ~ 1^- The form factors ^{chip) of 

the _7"^ = 0+, 1"*" baryons with the heavy quark spin switched off are called 
Isgur-Wise Q functions. For the form factors of the physical spin symmetry 
multiplets Aq, Eg, Eg we have from ( pT| ) 

{Bi\J\B2) =uiT[TT'2U2ac^ip). (20) 

This is equivalent to the results of 0,^. The contribution to the double 
spectral density is ^fif2^{chip)S{u)i — £i)5(w2 — £2)- 
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We adopt the standard continuum modeh the continuum spectral den- 
sity is equal to the theoretical one ([T^-[ig|) starting from a smooth curve — 
continuum threshold. Note that a non-smooth behavior of the continuum 
threshold imposed in the first paper of Jisf has no physical justification, and 
leads to an infinite slope of form factors at the origin. We choose the simplest 
variant — a straight line continuum threshold (Fig. This triangular con- 
tinuum model works well in the pion form factor case p^ . If the threshold 
is curved, then the contribution of the shaded regions in Fig. ^ should be 
subtracted, this will infiuencc the form factor slope. This degree of freedom 
for the slope is analogous to the freedom of varying the continuum threshold 
in two-point sum rules, and should not be very significant if the sum rules 
are applicable. With the straight line threshold at ti = T2, we don't need 
p{u>i, L1J2) to write down the sum rule; the one- variable function p{uj) (given 
by ( IitI) without the limit sign) is enough. Moreover, this function is propor- 
tional to where n is evident from the dimensional analysis. Hence spectral 
densitie s (p!8|-[19| ) are unnecessary: we can use terms of the coordinate space 
results (|13|-|14|) at ti = t2 = — 2r/2 multiplied by the corresponding fn{iL>cT). 

Using the dimensionless variables, we obtain the sum rules 



n^(ch ip)e 



-E,./E 



MEc/E) 



E% h{EjE) 
E^ ch^ ^ 



-exp(-2§ch^| 




f2iEc/E) 



§fo{E./E) . 



(21) 



agdch ^ 



27ttE^ 

Here c = 1 for A and c = — i for E; in the diagonal sum rule S = 1 for A2, 
S = —1 for and S — otherwise; in the nondiagonal sum rule d = chip for 
S|| and d = 1 otherwise. In the diagonal sum rule we have included the effect 
of noncoUinearity in the nonlocal quark condensate to the exponent; this is 
an order-of-magnitude estimate only, and should not be used when this effect 
is significant. At 1^9 = these sum rules coincide with ([l0|). Dividing ( ^ ) 
by ( |lO| ) , we obtain the formulae for the Isgur-Wise form factors not containing 
Er] the normalization at = is automatically correct. 
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Figure 7: Sum rules for the Isgur-Wise form factors, a) Ag: lower curve — 
the diagonal sum rule 1 at Ec = 4.6, E = 1.2; middle curve — the diagonal 
sum rule 2 at the same values; upper curve — the nondiagonal sum rule at 
-Ec = 3.9, £^ = 3. b)SQ'': lower curve the diagonal sum rule 1 for E|| at 
Ec = 5.6, E = 1.2; middle curve — the diagonal sum rule 2 for S|| and both 
diagonal sum rules for E_l at the same values; upper curve — the nondiagonal 
sum rule for S|| at Ec = 5.6, E = 3 (for T,± it differs by less than the line 
width) . 
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The results of the sum rules analyses for ^(chi^) are shown in Fig. 0. We 
have taken the best values of the continuum threshold Ec and the intervals 
of the Borel parameter E from the two-point sum rules. Variation of these 
parameters in reasonable bounds leads only to slight changes of these curves. 
The diagonal sum rules are very stable with respect to varying E, but only 
at -B > 1.2; this bound is higher than in the two-point case. This can be 
explained by a not very accurate consideration of the noncoUinearity effect 
in the nonlocal quark condensate. In the case of Ag, the two diagonal sum 
rules are in a reasonable agreement with each other; the nondiagonal sum 
rule predicts somewhat lower slope and a more straight shape, but is not 
in a sharp disagreement with the diagonal ones. In the case of Sg , the 
second diagonal sum rule predicts the equal form factors for the two possible 
light fields' polarizations; the first one yields the same result for S_l, but a 
somewhat different result for Therefore the accuracy of the approach 
doesn't allow us to distinguish between ^||(ch(y9) and £^±{chip). The results of 
the diagonal sum rules are in a reasonable agreement with each other. 

On the other hand, the nondiagonal sum rule predicts a significantly lower 
slope and a more straight shape, the nondiagonal sum rules for E|| and T.±_ 
differ only in the very small d = 9 (qq)^ term, and the predicted curves 
are indistinguishable. We remind that in the two-point case an agreement 
between the nondiagonal sum rule and the diagonal ones for Sg also was 
much worse than for Ag . It would be interesting to understand the reason of 
this poor agreement. 

Note that in the case of the nucleon form factors a thorough analysis of 
the sum rules appeared impossible, and a simplified local duality approach 
was used | |2^ . It corresponds to working at infinite Borel parameters, and 
using the continuum threshold from the two-point sum rules. Moreover, only 
one correlator was used, so no self-consistency check was possible. Here we 
are in a somewhat better position: we do have wide stability regions in Borel 
parameters, and comparison of different correlators allows us to estimate the 
accuracy (though it is not high). 
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